Scalar and Vector Perturbations in Quantum Cosmological Backgrounds 
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Generalizing a previous work concerning cosmological linear tensor perturbations, we show that 
the lagrangians and hamiltonians of cosmological linear scalar and vector perturbations can be put 
in simple form through the implementation of canonical transformations and redefinitions of the 
lapse function, without ever using the background classical equations of motion. In particular, if 
the matter content of the Universe is a perfect fluid, the hamiltonian of scalar perturbations can be 
reduced, as usual, to a hamiltonian of a scalar field with variable mass depending on background 
functions, independently of the fact that these functions satisfy the background Einstein's classical 
equations. These simple lagrangians and hamiltonians can then be used in situations where the 
background metric is also quantized, hence providing a substantial simplification over the direct 
approach originally developed by Halliwell and Hawking. 

PACS numbers: 98.80.Cq, 04.60.Ds 



I. INTRODUCTION 



In the theory of linear cosmological perturbations, sim- 
ple evolution equations for the perturbations have been 
obtained p). Lagrangians and hamiltonians describing 
the dynamics of scalar, vector, and tensor perturbations 
coming from the Einstein-Hilbert lagrangian have been 
greatly simplified in different cosmological scenarios un- 
der the assumption that the background metric satisfies 
Einstein classical field equations, and after taking out 
space and time total derivatives pj. Once these simple 
lagrangians and hamiltonians are obtained, the quanti- 
zation of linear cosmological perturbations becomes easy, 
with a quite simple interpretation: they can be seen as 
quantum fields which behave essentially as scalar fields 
with a time dependent effective mass. The time vary- 
ing background scale factor which is responsible for this 
"mass" acts as a pump field Q, creating or destroying 
modes of the perturbations. In this framework, one can 
assume an initial vacuum state for the perturbations, 
yielding primordial perturbation spectra which can be 
compared with observations. In the cosmological infla- 
tionary scenario |3| , the resulting spectrum for scalar per- 
turbations is in good agreement with the data [^|. 

However, this state of affairs is rather incomplete: the 
overwhelming majority of classical backgrounds possess 
an initial singularity at which the classical theory is 
expected to break down, and one needs to justify the 
initial conditions for inflation and quantum perturba- 
tions. Hence, a full quantum treatment including the 
background must be constructed. The first approach in 
this direction was made in Ref. where the canon- 
ical quantization of the perturbations and background 
was implemented through the derivation of the super- 
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hamiltonian constraint of the whole system and its conse- 
quent Wheeler-DeWitt equation H(A, P A ,X, P x )^ = 0, 
where A and Pa represent the phase space background 
variables, and X and Px the perturbation phase space 
variables. They claim that the no boundary proposal 
can set the initial conditions for inflation and the vac- 
uum initial state for the perturbations. Then, through 
the imposition of the ansatz on the wave functional 
*&(A,X,t) = (p(A, t)ip(A, X, t), they could manage to 
separate the quantum effects in the background from the 
quantum perturbations, where the wave function for the 
background tp(A, t) obeys an independent quantum min- 
isuperspace description where back reactions terms from 
the quantum perturbations are negligible. The singu- 
larity is bypassed through an euclidianization of space- 
time near it, and a consequent beginning of time when 
(or where) the geometry passes from the euclidian signa- 
ture to the lorentzian one. The quantum perturbations 
are described in the oscillatory part of the background 
wave function, where a WKB approximation can be used. 
Then, the evolution of the scale factor in time may be ob- 
tained through the equation a oc dS/da, where S is a so- 
lution of the classical Hamilton-Jacobi equation. Hence, 
this evolution is the classical one, and we are back to a 
semiclassical description of the perturbations. 

In parallel to that, the possibility that the singular- 
ity could be avoided through a bounce connecting the 
present expanding phase with a preceding contracting 
phase has been explored. In this case, the Universe 
is eternal: there is no beginning of time, nor horizons. 
Many frameworks where bounces may occur have been 
proposed 0, Q, @ ■ These new features of the background 
introduce a new picture for the evolution of cosmologi- 
cal perturbations: vacuum initial conditions may now be 
imposed when the Universe was very big and almost flat, 
and effects due to the contracting and bouncing phases, 
which are not present in models with a beginning of time, 
may change the subsequent evolution of perturbations in 
the expanding phase. Because of that, the evolution of 
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cosmological perturbations in bouncing models has been 
cause of intense debate |{| • 

In the framework of quantum cosmology in minisuper- 
space models, bouncing models had also been proposed 
where the bounce occurs due to quantum effects in the 
background [H EJ H 111 13 111 Hi Some approaches 
have used an ontological interpretation of quantum me- 
chanics, the Bohm-de Broglie one, to interpret the 
results HH)- In this interpretation, quantum 
Bohmian trajectories, the quantum evolution of the scale 
factor a q (t) at zeroth order, can be defined through the 
relation a cx dS/da, where <S* is now the phase of the back- 
ground wave function ip(A, t) without any approxima- 
tion: it is not a solution of the classical Hamilton- Jacobi 
equation. In fact it satisfies a modified Hamilton- Jacobi 
equation derived from the Wheeler-DeWitt equation for 
ip(A,t), and hence a q (t) is not the classical trajectory: 
in the regions where the quantum effects cannot be ne- 
glected, the quantum trajectory a q (t) performs a bounce 
which connect two asymptotic classical regions where the 
quantum effects are negligible. One than has in hands a 
definite function of time for the background, even at the 
quantum level, which realizes a soft transition from the 
contracting phase to the expanding one. Due to the re- 
sults of Ref. where the background minisuperspace 
Wheeler-DeWitt equation for <p(A, t) continue to hold 
when quantum perturbations are present because back 
reaction terms are negligible (which can also be justified 
through other ansatz for the wave function or verified 'a 
posteriori'), this background quantum function a q {i) is 
sufficient to describe the whole quantum features of the 
background. The natural question to ask is what hap- 
pens with the perturbations when it passes through this 
well defined and regular quantum bounce. One could 
then use the hamiltonian H of Ref. Q to investigate 
the evolution of quantum perturbations in this quantum 
background. However, the structure of H is rather com- 
plicated, turning it difficult to obtain any detailed result 
about the spectra of perturbations, specially the scalar 
ones. Also, a simplification of H using the zeroth order 
classical equations, as done in Ref. and described in 
the beginning of this section, is not possible because the 
background is also quantized and it does not satisfy the 
classical Einstein's equations. This state of affairs mo- 
tivated us to find a way to simplify the hamiltonian of 
Ref. p| , without ever recurring to the background classi- 
cal equations, and apply it to these quantum systems. 

Recently, we have managed to put the hamiltonian of 
tensor perturbations into a very simple form through the 
implementation of canonical transformations and redef- 
initions of the lapse function only, without recurring to 
any classical equations of motion |18| . Its consequences 
were explored in Ref. 0|. However, tensor perturbations 
are very special (they are automatically gauge invariant, 
their equations do not depend on the matter background) 
and it remained to investigate if it would be possible to 
do the same procedures to simplify the hamiltonian of 
scalar and vector perturbations. Note that such pertur- 



bations are not gauge invariant from the beginning, and 
they have contributions from the matter perturbations, 
which renders the calculations much more intricate. 

The aim of this paper is to show that it is indeed pos- 
sible to put the complicated hamiltonians of scalar and 
vector perturbations of Ref. .5] into the very simple form 
of Ref. without using any classical background equa- 
tions. We will exhibit the canonical transformations and 
lapse functions redefinitions which make the job. The 
simplified constraints obtained have direct physical in- 
terpretations. The quantization of the theory yields a 
very simple Wheeler-DeWitt equation for the perturba- 
tions and background, which can be used with whatever 
interpretation and choice of time one makes. 

When the matter content is a perfect fluid, the 
Wheeler-DeWitt equation assumes a Schrodinger form, 
and a further simplification can be achieved provided 
one uses the ontological interpretation of Bohm and de 
Broglie E3. As in this case a quantum Bohmian tra- 
jectory a q (t)) at zeroth order can be defined, a time de- 
pendent unitary transformation can be implemented in 
the scalar perturbation sector using this a q (t), and, like in 
Ref. , the hamiltonian can be further simplified render- 
ing equations governing the scalar perturbations which 
are formally equivalent to simple equations for a scalar 
field with an effective mass depending on the quantum 
solution for the scale factor a q (t), the quantum "pump 
field". 

The paper is organized as follows. In the following sec- 
tion, we specify the action and hamiltonian by restricting 
attention to the particular case of a Friedmann-Lemaitrc- 
Robertson- Walker (FLRW) background and perturba- 
tions around it, without yet making any separation in 
scalar, vector and tensor perturbations. In sections III 
and IV we analyze the cases of vector and scalar per- 
turbations, respectively. We concentrate in the hydro- 
dynamical fluid case letting the scalar field case to Ap- 
pendix A. In section IV we present all the steps to sim- 
plify the scalar part of the hamiltonian. In section V we 
quantize this system. After separating the background 
Schrodinger equation from the perturbed one, we show 
how to use the Bohm-dc Broglie interpretation in order to 
perform the last canonical transformations which yields 
quantum equations for the perturbations with the same 
form as those presented in Ref. |lj . Finally, Sec. VI ends 
this paper with some general conclusions. Appendix B 
presents the explicit canonical transformations used in 
section IV. 



II. LINEAR COSMOLOGICAL 
PERTURBATIONS 

Let the geometry of spacetime be given by 

9^ = 9$ +h llv , (1) 
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where g$ represents a homogeneous and isotropic cos- 
mological background, 

ds 2 = g^dafdx" = N 2 {t)dt 2 - a 2 {t) llJ dx t dx j , (2) 



Substituting Eq. Q into the Einstein- Hilbert action 
1 
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where 7^ is the spatial metric of the spacelike hypersur- 

faces with constant curvature K = 0,±1, and the h„ u where P = 8 ^ G / 3 > y ields the zeroth and second order 
are the linear perturbations, which we decompose into 
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The first order action was discarded because we are as- 
suming that the mean value of the perturbations over the 
spatial sections are null: 



d^xj^h^ = 0. 



(7) 



These are the actions for the gravitational sector. Let 
us now focus on the action for the matter sector. We will 
concentrate on perfect fluids for two reasons: 1) When 
quantizing the theory, a time variable appears naturally 
putting the Wheeler-DeWitt equation into a Schrodinger 
form. 2) The background quantum solutions of these 
models are quite simple [l3l Il5|. contrary to the scalar 
field case, where the quantum trajectories are implicity 
functions of time 01- We will let the discussion con- 
cerning the scalar field to Appendix 1, which follows the 
same lines as below. 

We will restrict the construction of the hamiltonian to 
the K = case and postpone the K ^ 0, which is more 
intricate, to a future publication. 

Following the approach of Ref . , the lagrangian den- 
sity of the perfect fluid is 



where 



and 



p — 



p[m + n(p, p)], 



p dp dp 1 p 
n dp' p' p 



(8) 
(9) 
(10) 



where p is the number density of particles, mo is their rest 
mass, and p is an arbitrary function of p, which will be 



identified with the pressure. The particle number density 
p is given by 



(ii) 



where F is an arbitrary function of lagrangian variables, 
a is a time parameter along the particle world lines, and 
J is the jacobian of the transformation from lagrangian 
variable to eulerian ones. 

The energy-momentum tensor of the fluid reads 

r" = ^ S (V^9 £ ) = £ yp.yv _ p ( g i*v _ V ^V U ), 
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( 12 ) 

where it is clear that e e p corresponds to the energy 
density and pressure, respectively. The sound velocity c s 
is defined by 



de 
dp 



(13) 



Perturbations displace the particles from their back- 
ground positions x^ to the x^ positions given by 



r(* ), 



(14) 



meaning a change into their eulerian position, which im- 
plies modifications in the jacobian, 



in the determinant, 
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(15) 
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The particle number density at point is then given by 
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Substituting all that in Eq. lfTU|) and finally in Eq. @, yields 



1 



drxNa?^ 

2 



1 , 1 



^ 2 + Ia'A: - 4>e \i ) ♦ -«> + r'-'< 



i 



2 -r ■ 4 - -.J 
1 



2 v- u , x-^y^s *. . N - 

The total lagrangian including the gravitational sector then reads 
d 2 aV 
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The procedure of Ref. |l| to simplify Eq. (121 )[l begins as follows: using the background equation of motion 

d 2 diV 3Z 2 7V 2 a 



2a N 



-Po, 



(21) 



and discarding a total time derivative 
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Now using the other background equation 



we arrive at the simplified lagrangian 
d 2 aV 
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Note that it is not necessary to use Eq. i|24[) in order to pass from Eq. 123(1 to Eq. 1(25(1 : the redefinition of the lapse 
function 



N =: N 



1 + w I rf3x7 " U! 



AM, 



(26) 



takes Eq. ((23(1 into Eq. I(25|) . Note that these two lapse functions related by Eq. ((26|) are equivalent at first order. 
Hence, this procedure does not modify the equations of motion at first order wnen we make a time gauge choice. 
Let us now calculate the hamiltonians of these lagrangians for perfect fluids with equation of state 



Po = Ae 



(27) 



The hamiltonian from Eq. 1(20(1 reads 
Nl 2 P 2 „P T Nl 2 P 2 
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N 
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while that from Eq. (|23|l is given by 



+ Na 3 e / d 



i xl h Lltf + I AM, - <PC |^ ) + Na 3 Xe J d 3 x^ Qe0 - 0f 



+ ^-Na 3 X(X + l) Eo / d 3 x 7 i [ -e 2 + C \i? \j ~ ,< ) (29) 



The quantity Pt appearing in the second term of the 
zeroth order term of both hamiltonians is just the kine- 
matical constant Pt = £oa 3A+3 U. We have introduced it 
as a canonical momentum to a variable T which is cyclic, 
implying indeed that Pt is a constant. We have made an 
inverse Routh procedure. The variable T plays the role 
of time when the system is quantized. This form of the 
zeroth order hamltonian appears in other approaches to 
a lagrangian formulation of fluids; see e.g. Ref. [2(J for 
details. 

One can now use the total time derivative (1221) to con- 
struct the generator of canonical transformations 



T = aP a 
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Using the fact that p oc a™ 3 , the particle number density 
transforms to 



P = P 



1 - 
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Substituting this last equation into Eqs. (|1U[) and JHJ we 
obtain 



(e~o + Po) x 
£0 = £0 ; op. 



from Eq. to Eq. I|25|) just through a redefinition of 
N, then it is proven that lagrangian (|20|l is equivalent to 
lagrangian l|25|l at this order of approximation irrespec- 
tive of the classical background equations of motion. 

In order to proceed from this point 1 , we will now sep- 
arate the perturbations into scalar, vector, and tensor 
perturbations. We make the decomposition: 

A, = B {i + Si 

dj = 2ipjij - 2E\ i \ j - F il:j - + wu (34) 

in the gravitational sector, while the quantities Wy, P; e 
Si satisfy 



S" 



F l m = 



and 



with 



«,« y = 

w\ = 0, 



V* \i = o, 



(35) 



(36) 



(37) 



(33) 



in the matter sector. Substituting the above decomposi- 
tions into eq. 1251) leads to a separation of this lagrangian, 
with some total derivatives discarded, into three indepen- 
dent sectors: scalar, vector and tensor sectors. We will 
focus our attention on the vector and scalar sectors be- 
cause the case of tensor perturbations has already been 
treated in Ref. [is| . 



Inserting Eqs. (0TJ and O into Eq. J5SJ), we obtain 
Hence, in the lagrangian point of view, one can pass from 
Eq. ij2TH) to Eq. (|2*3*|> without using any background equa- 
tions of motion. As we have shown that we can pass 



1 Equation 1251 corresponds to Eq. (10.37) of Ref. y| if one is 
restricted to scalar perturbations, and if one reads j3 in the latter 
as /3 = 3a 2 i 2 (e +po)/2. 
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III. VECTOR PERTURBATIONS 



and 



Combining the contributions of gravitational and mat- 
ter sectors and defining the gauge invariant quantities 



rj l M = rf + F\ 



(39) 



V 1 = S* 



N 



-F l 



(38) we obtain 
I 



(40) 



When constructing the hamiltonian, the primary constraint II, « appears, where ss means a weak equality in the 
sense of Dirac |2l|, and Ilj is the momentum canonically conjugate to V 1 . The hamiltonian then reads 



H v = N 



IP 7, Pi 



iaV a 3X 



d A x 



P l P, 



V l P, 



p Ta 2-3A 7 i/2 a 121 



a i 

7 2 



yi\iy. 



>\j 



d 3 xAV, : 



(41) 



where Pi is the momentum canonically conjugate to if ( gt ', and the A 1 are Lagrange multipliers. 



The conservation of H in time imposes a secondary 
constraint 



IF = {U\H V } = 4>\ = -tt* - b , (42) 



The two constraints obtained are second class. After 
defining the corresponding Dirac brackets |2l|, they be- 
come strong equalities which can be used to obtain some 
variables from others. 



The conservation in time of fixes the lagrange multi- 
pliers A 1 to the value 



l 2 P 
A 1 = —2-V 1 
a 2 V 



(43) 



Then the equations of motion for V and rf ^ imply 
that 



IV. SCALAR PERTURBATIONS 



Defining the quantities 



F = B E, 

N ' 



(46) 



V* = -f , (44) 

and 

,< s (-. ( ,, + ^ = _v!|_ (45) 

These solutions correspond to the classical result, which 
was obtained without recurring to the classical back- 
ground equations. 



C G?i) = £ + E> (47) 
which is a gauge invariant quantity, and 

y= ^aVjA + l) g0 ( a- (g , +n (4g) 

which can be identified with the perturbed velocity po- 
tential of the fluid particles, the scalar lagrangian reads 
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d 3 xj 2 (ij; + 
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(49) 
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As in the vector sector, some constraints appear, and, because of definition (|48|l which involves a time derivative, 
we have to use the Ostrogradsky method |22( through the definition ir^ = dLjdC, — tt^. The constraints are 



k = Pn; fa = ?tf; fa = fa = n v ; fa = P M , 



(50) 



and the hamiltonian is 



H = N 



H a + J d 3 xA 4> ir <t , + J d 3 xA F ir F + J d 3 xA lp n if 



AnP, 
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where Hn reads 
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(52) 



Conservation in time of the primary constraints 150(1 
leads to the secondary constraints 



H 

N 

1 2a if I 2 

(X + 1)P T i 



a 3X V 

i Vv 



The Poisson brackets among the constraints read 
B (sO} = -^7*<^-^« 

{0 7 (x), <?!>8 (s')} = -TH~l^ 5 i x ~ X 'Y ,i 

2 i I 2 

{fa(x),fa(x')} = —^8{X-X'Y J tij - —TT i ,(x)TT <: (x') 

{fa(x),fa(x>)} = ^J-^ a-i^Mfo - . 



6/ 2 a 



(53) 



4V6v/y(A + l)P T 



(1 - 3X)a-^ 1+3X K 4 ,(x)Tr <: (x'). 



(57) 



Neglecting third order terms, conservations in time of fa 
and fa are identically satisfied, whereas conservation of 
fa, and fa determines the lagrange multipliers Ap and 
Am. The Lagrange multiplier Ap reads 2 



Defining 



»6 T ~<P2 

a 



l A , V(A + l)P T ^ a _i (1+3A) 

VaVT 7 



07, (58) 



l 2 P 

aAp — tjj — cf> -\ 7j-F- 

av 



(54) 



As A F = F/N, then 



N 



F = ip ■ 



2d 



-F. 



which, when expressed in terms of the gauge invariant 
Bardeen potentials, yields 



$ = * 



a well known result. 



one can prove that fa is a first class constraint: it has 
zero Poisson brackets with all others constraints up to 
third order. We are then left with four second class con- 
straints, fa, 05, 07 and 08, Hence, from the 10 degrees 
of freedom of phase space corresponding to the variables 
(55) 0, ip, F, (p and £, we have to extract 4 from the second 
class constraints and 2x2 = 4 from the two first class 
constraints fa and fa , remaining 2 degrees of freedom in 
phase space, as expected for this problem. 

In order to eliminate the second class constraints, we 
have to define the Dirac brackets associated with them. 
The Dirac brackets among the variables of phase space 
which are not canonical are (excepting the ones involving 
F and ttf, which are not relevant) 



(56) 



The explicit value of K v is not important for what follows. 



{C^f{x)^{x')} D 



VeiWa-^ 1 - 3 ^ 



W(A + 1)Pt7 



-S(x — x') 
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{P a tp(x)} D = —(l — 3X)ip(x) (59) we obtain that the Dirac brackets for these quantities 

2a are canonical. The hamiltonian in terms of these new 



Defining the quantities 

nc) =: a ^- 3A V, 



variables then reads 



vV (A + l)P T . . Wl H = NH„-Nj d'xM, + J tfxK^ (61) 



3-\/A\/(A + 1)P T i , . 

i"V(c) = 7i> / -, t= 7 2 V(c), (60) where i? is given by 

Votv v 

I 



l 2 P 2 Pt 3P 
4aV a 3A a 3A 

VA /3V(A + l)iV 



2 V 2 a 3Vy 
and 06 by (we omitted the bars) 



A 



12Z 2 a 2 ~ 



3A(A + 1)P T 
8a 3 y 



d 3 x 7 ^ 2 c) 



d xtp^nj, (c) + I d ,i 



V6V(a + i)Pt _ 3A ; 2 P a 

-a 7T V ( c ) 



7Ti/> (c 



VICJ 2a 2 ^' 

From the second class constraints we obtain the identity, 



3VALP aV /(A+l)P T 
2 v / 6a 2 ^t 



i 2a i 



2^^ + ^ ■* = 
which in terms of the new canonical variables reads 

I 2 



3l 2 VX^/{\ + 1)P T 
2 y/6ly/V 



2a 2 72 



tt«/> (c) + g-F 1 ' 1 ,i = 0. 



We will need this equation later. 

If we now perform a canonical transformation generated by 



Ft = TP T + <iP„ / d 3 x 

where a is given by 



1 



a 2 (1 3A) ^( C )7T + Vwv> + 



2Vy v /(A + i)p T a 



l 2 P a VX 



(A + l)Pr (l-3A)P n „_ (2 _ 
2Z 2 P a a 24V 



(2-3A) 



(62) 



(63) 



(64) 



(65) 



(66) 



(67) 



constructed in order to introduce the v variable of Ref. Q] (ir is its canonical momentum) 3 , the new Hq reads (the 
explicit canonical transformations are given in Appendix B) 



H 



VP 2 P T 1 



AaV a 3X 2a 



A 



+ — d 6 x— + — d a x^v' % v ti + 



7 . 



2a 



9[(A + 1)P T ] 2 (1+6A) _ 9A[(A + l)P T ]P T (1+6A) 
2P 2 " 2P 2 



(-4+18A-18A 2 )Z 4 P 2 6Z 2 (1 + 3A 2 )P T 



64a 3 ^ 2 



-(2+3A) 



16V 



d 3 x7 2 w 2 + 



^ 2 Pa 



2 (1 3A) / (f x ^ 2v ^,: 



3 The term proportional to a in Eq. 1661 . as well as the specific 
form given by Eq. 1671 are made in order to eliminate a term 
proportional to vtt in the hamiltonian. 
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V 



18[(A + 1)P T ] 3 ^ (2 _ 9A) _ 18[(X + 1)P T ] 2 P T V „ (2 - 9X) 



(_ 2 + 9A-9A 2 )[(A + l)P T ] _ 3A 3[(A + 1)P t ] 2 (1 _ 6A) 
8A^ Ppa A 

(-2 + 3A) v /(A+l)PT a _i a+3A , _ 6A^V(A + l)P r PT a x a _ 9A1 
2VA\/V F / 2 p„ 2 Va 



^4P4A /4p4 

3(3 + 2A + 3A 2 )[(A + l)P T ]Pr n 



(1-6A) 



2l 2 P 2 X 



72 ^ 



+ 



18[(A + l)P T ]^ a i (1 _ 15A) 18A[(A + l)Pr]'PrVF a i (1 _ 15A) (2 - 9A + 9A 2 )Z 2 P aV /(A + l)P T a _3 (1+A) 
/ 2 P Q 3 VA ' PP a 3 VA " 8V*VA 



3[(A + 1)P T ] 



i a -i(i+9A) _ 3(3 + 2A + 3A 2 ) ^(A + 1)P T P T a _ |(1+9A) 



VXy/VP a 

_a_ _ 2(A+l)P T y 2 _ 3A 

3/ 2 ; 4 p 2 



1 , , 3(A + 1)P T _ a+3A) 

a 



7 2 v 



(68) 



This canonical transformation applied to (f>e shows that v is a gauge invariant quantity. The same is not true for 
its momentum tt, which has a non-zero Poisson Bracket with the first class constraint <p e . In order to obtain a gauge 
invariant momentum tt we now make the canonical transformation generated by 



T2 = aP a + J (fx{4>iT^ + VTT + 

6F[(A+1)P T ] 2 3(1 _ 2A) (1 + 3A)(A + 1)P T 2 _ 3X 



6y^[(A + l)Pr]i a 3 fl _3 A) _ (l + 3A)V(A + l)P T aA(1 _ 3A) 
l 2 PaVX 



; 4 p, 3 a 



2/ 2 P„A 



1 ,2 2a3V , , i 

7 2 ^ + -TTTri 3 ^ ,i 



H 4 Pa 



(69) 



and aiming at eliminating a term in v 2 proportional to 
Pt in the final form of the hamiltonian, we perform the 
last canonical transformation 

Tz = aP a + - [ d^xvir - [ d 3 x^v 2 . (70) 

a J 4aV J 

The constraint </>6 reads then 

l 2 P a 



and as — is not weakly zero, we can redefine the 
constraint </> 6 to be 

06 = 7TV ( 72 ) 

This constraint, in the Dirac quantization scheme, will 
imply that the wave functional will not depend on ip, 
and the second class constraint ijB^ll turns out to be the 
usual relation between cf> and v from Ref. pj, Eq. (12.8), 
as we will see later on. The new Hq is given by 



)2p2 



#0 



1 



AaV a 3X 2a 3 



PI , Pt 



,0 tt aX 
<Fx— + — 



d xj 2 v' l v,,, 



^ 3lV(A + l)f T | (1+A1 3(A + 1)P T (1 _ 



■3A), 



H, 



(0) 



d 3 x< - 



V 

9(X+1)P T ! 



P, 



ip> d XTT^ 



a 1 - 3A 7 ^ 2 + 



VPl 

18[(A + 1)Pt]*-v/V 



6%/^ y/(A + 1)P T 



pp 2 ^/x 



PP2VX 



2a 2 V 
l^Pj 



ijjip' 



„f (1-8A) , 9(A-l)y/(A + l)Pr i a _ 3A) 
2VF\/AP a 
9(l + 2A)[(A + l)Pr , 3A 9A(A-1)^ 2 



3X ^TT 

72 <mj) 



H, 



(o) 



(2) 



d 6 x 



2P a 2 

I 



81" 



(73) 



where Hq ^ and are the zeroth and second order second order functions which can be read from Eq. I|73|l . 

hamiltonian constraints, and F^ 1 ) and are first and 
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We now make the redefinitions of N and cj) as N = 
N(l + J d 3 xF^), which would again just imply a differ- 
ent irrelevant time gauge choice with terms beyond first 
order, and 4> — N(— ^Sv ^ + F^). From the inverse of 



and 



the transformations (|70(l and (|60|l . and definition (|48|l . 
is given by 



l 2 P 
2a 2 V' 



3/ 2 n ■ 

—Na(\ + l)e (-t + B) (74) 



As <j> is, through the equations of motion, equal to ip, we 
obtain, imposing N = a, the constraint equation (10.39) 
of reference 0. 

Inserting expression l|73[l into Eq. (|61[) . and the above 
redefinitions, we obtain, omitting the tilda, 

H = N(H { 0) +H^ 2) ) + A N P N + J d 3 x<jxt> & + J d 3 xA <t ,ir 4> , 

(75) 

(76) 



with 



(o) = 



AaV 



Pt 

„3A' 



1 

2^3 



d 3 3 



7 5 



a A 

T 



(77) 



Now we are left with two first class constraints (in 
fact one plus oo 3 constraints): one with the homoge- 
neous lapse function N as its associated Lagrange multi- 
plier, which in the quantization procedure will lead to the 
Wheeler-DeWitt equation, and the other oo 3 constraints 
with 4>{x l ) as their Lagrange multiplier, which is noth- 
ing but the inhomogeneous lapse function (see definition 
(|3*|)). which, as anticipated, has been tremendously sim- 
plified to imply a simple consequence when quantized: 
the wave functional does not depend on ip. The super- 
momentum constraint is automatically satisfied because 
the v variable is gauge invariant. 

The connection between v — av (the Mukhanov-Sasaki 
variable) and $ can be obtained from Eq. I|65(l which, af- 
ter implementing the canonical transformations l|t)6l I69|) 
reads (the bars are omitted) 



2a^ ' 3l 2 P a 

+ #TWy fl4(1+3A) ^ o 



2 " V ^ + ^F 
V 2aV 



3W{(\ + l)P T }? a _i (1+9X) _ (l + 3A) v /(A + l)F T / 2 a _ 3/2(1+A) 
P 2 V\ ' 4y/\VV 



(78) 



Using that P a 
i 

3Z 4 P a 



2Vad 



N[2 , we can identify the quantity ip + ^fiF with the Bardeen potential ^ which, from Eq. i|56[l . 
is equal to <I>. After some algebraic manipulations, we obtain: 



r 7ty + & 



3/V(A + l)P T a _3 (1+A) f 7T 



l 2 P a , 3(A+1) 



2aV 



P, 



H )v 



0(3) = 



(79) 



Using again that P a = — 2 ^ a , and that n = 72 v, Ho « 
0, 7T.0 0, and choosing the gauge N = a (conformal 
time), we obtain 



3Zy(A + l)Pr a _x (1+3A ^Y 



2VaVu 



(80) 



Equation ijSUJ) coincides with equation (12.8) of Ref. Q 
relating v and $ when the classical equations of motion 
are used. 



V. DIRAC QUANTIZATION 

In this section we will focus only in the quantization 
of scalar perturbations. Vector perturbations are trivial 



and the quantization of tensor perturbations was done in 
Ref. [3. 



A. The functional Schrodinger equation 

In the Dirac quantization procedure, the first 
class constraints must annihilate the wave functional 
X[N, a, <j)(x l ), %p(x),v(x l ), T], yielding 

d 

dN X = °' 
T7X = 0, 



H X = 0. 



(81) 
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The first three equations impose that the wave func- 
tional does not depend on N, (f> and i^: as mentioned 
above, N and <j> are, respectively, the homogeneous and 
inhomogeneous parts of the total lapse function, which 
are just lagrange multipliers of constraints, and tp has 
been substituted by v(x l ), the unique degree of freedom 
of scalar perturbations, as expected. 



As Pt appears linearly in H, and making the gauge 
choice N — a 3A , one can interpret the T variable as a 
time parameter. Hence, the equation 



H X = 
assumes the Schrodinger form 



(82) 



J 



dT X ~A\ da 



a (3A-l)/2ii 

da 



— / dx w ~ —7- ' dxvv < 



(83) 



r 



where we have chosen the factor ordering in a in order to where Xto){ a t T) satisfies the equation, 
yield a covariant Schrodinger equation under field redefi- 
nitions, and V and I have been absorbed in redefinitions 
of the fields. 



B. Further developments using the Bohm-de 
Broglie interpretation 



4 r da 



fl (3A-l)/2_^_ 

da 



X(0) (a,T), (85) 



If one makes the ansatz 

X[a, v, T) = x(o) (a, T)x(2) [v, T] 



(84) then we obtain for X(2) (o> v> T) the equation 
I 



(3A-1) r S 2 

/ d A x-— 

2 / 6v 2 



X(2)(a,v,T) 
I 



Aa (3A+l) 



d^xv^v^X(2){a,v,T) 



(86) 



Solutions of the zeroth order equation l|85(l are known 
[l3l fill . If one uses the ontological Bohm-de Broglie 
interpretation of quantum mechanics in order to obtain 
the bohmian trajectories a(T) from Eq. I|85(l . this a(T) 
can be viewed as a given function of time in the second 



equation l|85|) . Going to conformal time dr\ 
and performing the unitary transformation 



U 



,{«[/ d 3 xy? f|]} {*[/ d 3 x(=±IZ) i n( i) ]} 



aT, 



(87) 



the Schrodinger functional equation for the perturbations 
is transformed to 



. dX(2)[v,v] 



d 3 x 



IS 2 X 



2a 



The corresponding time evolution equation for the op- 
erator v in the Heisenberg picture is given by 



v )X(2)[v,v], 



where we have gone to the new quantum variable v = av, 
the Mukhanov-Sasaki variable defined in Ref. 0, after 
performing transformation l|87|) , and we have omitted the 
bars. 



Xv' 



0, 



(89) 



where a prime means derivative with respect to confor- 
mal time. In terms of the normal modes Vk, the above 
equation reads 



Vk+ [k -— )v k = I). 



(90) 



These equations have the same form as the equations 
for scalar perturbations obtained in Ref. (for one sin- 
gle fluid, the pump function z" / z obtained in Q is ex- 
actly equal to a" /a obtained here, if we make use of the 
background equations). The difference is that the func- 
tion a(rj) is no longer a classical solution of the back- 
ground equations but a quantum Bohmian trajectory of 
the quantized background, which may lead to different 
power spectra. 
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VI. CONCLUSION 

In this paper we have managed to obtain simple hamil- 
tonians for scalar perturbations when the matter content 
is described either by a perfect fluid or by a scalar fluid, 
without recurring to the background classical equations. 
Performing canonical transformations and redefining the 
homogeneous lapse functions with terms which do not al- 
ter the linear perturbation equations, the constraint con- 
nected to the inhomogeneous part of the lapse function is 
greatly simplified implying that the momentum canoni- 
cally conjugate to the scalar perturbation tp is weakly 
zero. The hamiltonian constraint is also greatly simpli- 
fied when written in terms of a new variable which is 
exactly equal to the usual Mukhanov-Sasaki's variable 

II • 

This simplified hamiltonian can now be used in 

the Dirac quantization procedure not only to quantize 
the perturbations but also the background, yielding a 
Wheeler-DeWitt equation much simpler to handle then 
the one of Ref. 0. In the case of perfect fluids, where a 
preferred time variable appears and the Wheeler-DeWitt 
equation can be put in a Schrodinger form, and using 
the Bohm-de Broglie interpretation of quantum mechan- 
ics to perform a last unitary transformation, one obtains 
an equation for the modes which has the same form as 
in Ref. , where the pump field is obtained from a scale 
factor which now takes into account the quantum effects, 
the quantum Bohmian trajectory of the background. 



In future publications, we will apply these results to 
specific models, and evaluate the power spectrum of 
scalar perturbations which arise on them in order to com- 
pare, when taken together with the results of Ref. [l9j . 
with observations. 



VII. ACKNOWLEDGMENTS 

We would like to thank CNPq of Brazil, CAPES 
(Brazil) and COFECUB (France) for financial support. 
We very gratefully acknowledge various enlightening con- 
versations with Patrick Peter. 

APPENDIX A: THE SCALAR FIELD 

In this appendix we implement the same simplifica- 
tions we have done for hydrodynamical matter for the 
case of a scalar field. 

The scalar field lagrangian reads 

£ m = \pwV* - \u(p). (Al) 
We write its perturbation as 

<p = ip + Sip, (A2) 

where <po is the homogeneous scalar field depending only 
on time. Substituting IA2I in IA1I we obtain 



£, 



2 A" 2 2 U 



N 2 ' 



2<A) 
N 2 



fa 
2N 2 



A l Ai 



"/V 2 " _ ~~Na |l 



1 



2N 2 



Sep 



■^S(p l S(pi - ^U vv S(p 2 



(A3) 



The total lagrangian including the gravitational sector reads 

a 2 aV NKaV tp 2 a 3 V Na 3 VU Na f j3 i r , ih -, 1 ll]k 

L = „ 1 s h— 1 77 / (Tarp lA^Aui,] e y| e i7 -i fc 

p N p ^ 2N 2 6Z 2 J ' 1 Nj| 4 lJlk 

+ ^ A ^ j tf + 2^ ^ k |fc + ^ I' " 2^ I' " ^ + 4 e|j£ ' 1 + K{ ~/ ~ ***** ~ 

_3 e 2 + 3 ^ + 3 eyey) _|a J d s xli{(j)A i ii _^ Aieijij) + ^_ J ^(6^-^6-06) 

~W j ^xj^tA* |,- — —j^- J d 3 x^(c/}+ i e )5ip + a 2 <po / d 3 xjhepA i ]i - Na ^n> J tf^iy _ I e )^ 

+ J d 3 ^ (30 2 + e4> - A* Al - ^ + \e 2 ) + J f^itf + e0 _ A'ij + I t «t« - \e 2 ) 

Na - 3 f j3 ±/Sip 2 Sip\ l Sip\i 1 2 

+ — J d x l 2 (^2 ^- L - (A4) 
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Its hamiltonian is given by 

l 2 P 2 KaV . Pi a 3 VU . l 2 P 2 f ± , . Z 2 P 2 /" ± Z 2 P 2 



H = Ni - 



AaV I 2 2a 3 V 2 8aF 2 

2p2 /• . ;2 p2 



/ - 3^ / ^ + W J ***** ^ + dSx ^ iJ w 



;2 /■ .. q;2 /• _2 ;2 



+-3 y d 3 *-^ _ _ y ^ _ __ y , 3x7re _ _ y ^ | . _ _ y ^ | 

/ 2 Pa /\ 3 , 2Z 2 P a /■ - .. 1 /■ - 7T 2 a /■ i 

y w y ^ eij + ^y d 3 ,-j - ^ y ^ 



+<t>\ie ij \J - \j - + \^S* + K{1 -/ - -z<P+ A ' A - - 

,37 



J d 3 x ^( < j ) 2 + ecf>-A i A i + ^e ij - 1 -e 2 ) + ^ j d 3 x^ (±6^6^ + \u w 6<p 2 ) + J d 3 xjhA\ 

(A5) 

Performing the canonical transformation generated by 



F = aP a + ip P v - J d 3 x(4>TT lP + A t ir\ + + 5ipw v ) - ^ J d 3 x^(^el J - ) 



^ J d 3 x^^+\e)5^ 
(A6) 



which are 



Pa = Pa-^f**rH&«J-\?) 

1 f 1 1 

ip = if + — d 3 xj^{4> + ^) 6( P 

6V 2 ' 2V 

^ = ^-^-yHt+l'e) (A7) 



yields the new hamiltonian 



l 2 P 2 a KaV Pi a 3 VU l 2 P 2 f j3 , 2 Z 2 P 2 f j3 , Z 2 P 2 , „ 

^ = ^ -737 - -J5- + 7777 + o + 7777 / ^^7^ 2 + 7^7 / d 3 X^e<f>- / d X-f^A z Ai 



AaV I 2 2a 3 V 2 8al/ 2 ./ ' r 8ay 2 ./ ' r 8a^ 2 

32 



3Z 2 P a 

+ 



^ y y ds^^,.--^ y d 3 X7 ^ 2+ ^_^ )+ ^ y d 3 OT0 

Jd^ + ^Jd^h^^ + ^Jd'x^-^Jd^-^JdWA^^^ 
-~a I J *M-^J ^, + ±J d 3 4 ,f.,„„, 

+a 3 U v J d 3 x 1 ^S^-^- J d 3 x^(0 2 +^-A tA t ) + Y J d 3 x^(^Sip\ l 5i PH + ^U vv S^ 2 )'j 



(A8) 
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Going back to its corresponding lagrangian, and redefining N as 

1 



N = N 



1 + — / d 3 xjHect> + <p-A x Ai) 



(A9) 



we obtain 



a z aV NKaV <p 2 a 3 V Na 3 VU Na f 3 i 



L 



l 2 N 



I 2 



2N 



+d>t-e tJ i et e lJ i ■ 



6Z 2 
1 



a 

N 



1 



+ + K(- -tPey + 2e<f> - 2A l Ai) 



24PN 



+ e \i e i \k 
3 



24Z 2 A 



d 3 x^i 2 



+ 



ad 2 
"PN 

aVo 

N 
Na 3 



d 3 xj2((j)A l 



A^ 



it \jj 



9 • 

a a 
3l 2 N 



d xj 2 (j)e 



, 3 i,^ 2 <M 1 <Vk l r , r 2 % 



6l 2 



d 3 x~f 2 eA l | j 



<A)V 
27V 



d 3 xj 2 4> 2 



(A10) 



Splitting as before the perturbations into their tenso- 
rial, vector and scalar parts 

Ai = B H + Si 



21/rfij - 2E W - F i{j - Fj\i + w l3 , (All) 



with 



» tJ M = 
= 



(A12) 



the lagrangian also splits in tensor, vector and scalar 
parts. The tensor part was already treated in Ref. |lSj . 



S l \i = F l K = 



The vector part reads 



L<" = dW>S^S m - ^ J d 3 x^S iU F^ + yj^ / dW-F^ + ^ / d»xr*^ 



6Z 2 iV / ' 3Z 2 



d 3 xjiF*F t - — — / (Far/* S 1 Si (A13) 



Using the gauge invariant quantity 

Vi = ^ - ^F,, (A14) 
this lagrangian simplify to 

^ | ^(l^V^ - 2AV*Vtf (A15) 
Its associated hamiltonian reads 

(A16) 

where we have the constraint 

ttV « (A17) 



Conservation of the constraint i|, « leads to the 
secondary constraint 

-y i|j '|j+is:Fi wQ (A18) 

whose conservation fixes the Lagrange multiplier Ay, 
which means that both constraints are second class. 
Defining the associated Dirac brackets, they become 
strong equalities, yielding the well known result for a 
universe filled only with a scalar field: 

V* = 0. (A19) 
In the scalar sector we have 
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Na 
2a 2 d 



^ + ^(-6^ + 12^) 



l 2 N 
Na 3 



,3 1/ 1 r - 2 1 . ir 1 . o, 2a 2 

(TxT(jp5<p --p5<p' % 5<p, i --U vip 5<p*)--y2 



^ a 2~ ^ 



2 a v N ' 



a 



(A20) 



When constructing the hamiltonian, we obtain the primary constraints (here again we define F — B — aE/N) 

<j>i=P N ~ 0, 



h = 7!> — 7 2 <V ~ 0, 



N 



TT F W 0. 



The hamiltonian reads 
where iJ is given by 

,2p2 

Ho = - 



H = NH a + A N P N 



J d 3 xA F ir F + j d 3 xA <j> {-K 4> 



P 



pi 



KaV a 3 VU l 2 P a , ,, 



p2 



AaV 2a 3 V I 2 
3 r .12 i 

2 fl 3 7 T^ + 3^'' < l)2 + ( ~ 2a 2 V 2 



_£_ J rf 3 ^. 

~W 2 j d ^ 

a f .o i 

- / d 3 XJ2 



2a 2 V J r " 2a 3 V 2 

3^ Pa Pip 



d 3 x^ 2 <f> + 



(A21) 



(A22) 



'■P / ^3 



2a 3 V 



d Sipir^ 



+ a A U v ) / d 3 x^(j)S(p 



- 2 </>Vi + ^(-3V^ 2 + 60V) 



'->/ 



+ 



Sip'Sipj + (- 



<»Z 2 ^ a 3 ^ 



4a 3 V 2 



) / d 3 x^5ip 2 



(A23) 



r 



Performing the canonical transformation 



where I represents the identity transformation, the new 
reads 



= V?o - J d 3 x^(j)8ip 



7T0 = + ^-^8(f 



generated by 



T = Z 



d xj 2 



(A24) 



H = NH Q + A N P N + J d 3 xA F TT F + J d 3 xA lP TT 4 „ (A26) 



(A25) where 
I 



Hn = - 



l 2 P? 



AaV 2aW 



p2 



KaV a 3 VU l 2 P a 



2a 2 V 



3/2p 



/ f I 1 



2a 3 7 5 
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+(- 



3l 2 P a P v a s U v 
2a 2 V 2 2 



) / cfarp 



a 



aK 



dTx^FF'- 



a 3 V 



d x6it w + 



2a 3 



3„K 



d 6 x 



+ 



9PPI 
4a 3 V 2 



d 3 xj 2 5ip' l Sip,i + (- 



+ 



a 3 £4 



) / d 3 xj 2 5ip 2 



(A27) 



Conservation of the primary constraints (|A21|I leads to 
the secondary constraints 



Pn«0 



(A28) 



2 i IK i 

{05, <M - - g727^ 3 (£ " x')-' 1 ,i ~ 3p7^ 3 (^ - x'). 

(A33) 



1 2a 

^3^ + 9^ ■< 



|f 7 ^F«0 (A29) 



; 2 P a ,3Z 2 P a P v a 3 [/^ i 



2a 2 ^ 
2a 



2a 2 V 
6aif 



3/ 2 



j 2 ip - 





(A30) 



Conservation of Po is identically satisfied. Conserva- 
tion of 06 leads to a term proportional to Hq up to second 
order terms. Finally 05 fixes the Lagrange multiplier Ap: 



a A, 



■0-0 + 



aV 



(A31) 



Substituting F — {F, H} — Ap into the above equation, 
we get for the gauge invariant Bardeen potentials $ and 



$ = * 



(A32) 



Calculating the non null Poisson brackets among the 
constraints yields 



The 0q and 05 constraints are second class, while 



"6 =• <P6 ~r ~<P3 
a 



(A34) 



is a first class constraint. Defining the Dirac brackets, 
the second class constraints can be substituted in the 
hamiltonian. 

Making K = 0, and performing the canonical transfor- 
mations generated by 



Fx = aP a + ip P v 



" - e- 2 



where 



d 3 o 



3P 2 



2P 

air Sip + tpTty - j2Tr' l P 7T 

' P a 

(A35) 



al 2 P, 



aP a V 2V ' 



(A36) 



{03,05} 



2a 
91 2 



j 2 5 3 (x - x'Y 



2aK i o , 
-^-7=<5 3 (x - a') 



and 



- 5 , / , 3 . , ,2P<^ a 4 VT^ 6P V i , 2P V a A VU v P v 6P V , i (2 
^ 2 = aP a + ^ P ip + / d 3 x 0>^+^+ — 7-— ^ ^_) 7 2^+ * J^JP ^3-7 2 ^ 2 

+7^7^'\i), (A37) 



3Z 4 P a 

and making a redefinition of AT, we finally obtain 

06 = Ttyj (A38) 

and 

/" / / 2 P 3P 2 3/ 2 P \ /" 

H = N Hq + J d 3 x f - ^0 + ^Ty^ + ^4T/ U ) & + AjvPv + y d 3 zA 07 r , (A39) 
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where 

l2P a P l a 3 VU 1 /\ 3 7T 2 1 [ - I t 

/15^ of^_3^a 9UPl_J^ 27P^_3P^\ f , 

+ \4a-V 2 + 4 8 4aP 2 16a 3 V 2 4 a 7y 2 P 2 P a J J 7 ' 1 j 

Using the background classical equations one can show that the coefficient of v 2 can be written as z"/z as in Q]. 
Without their use, this is the simplest form the hamiltonian of scalar perturbations can have in scalar field models. 

APPENDIX B: THE EXPLICIT CANONICAL TRANSFORMATIONS 

The explicit canonical transformations obtained from the generators T\, Ti and T-$ of section III are, respectively 



1+ V^t^^ d-aA) l^+jyto J d ^ X7 i(V6la^-^v- 

P a = P _(l^a-i(X-A) /^ ( ^ (1 -3^_ gVMA+WF 



2 v / 6V(A + l)P T yU. 2 _ 



6 2 - 3A ^) 2 



/^_ifl 3A 1 276^+1)^^.2 , A7 

/p Q vA 

2 v /(A + l)iWv r _3 (1 _ A) 
^ (c) = ^ + «>C V 



(Bl) 



, 12VF[(A + 1)Pt]t.i, 91 QX) f , 3 i - 
1 H !+= — _: J a? (91 ' M) / d 3 x75ui/' + 

pp„ 3 7a J 



18U[(A + 1)P T ] 2 A , (l-3A)(Al)Pr sl _3 A ' 



Z 4 P 4 A 



2a 2 V 



3Z 4 P 2 



Pa = Pa 



9(l-3A)7y[(A + l)P T ]i 5 i (1 . 9A) _ (l-9A 2 )7(A + l)P T a _: (1+3A) 

z 2 p 2 \/a ' 4\/A\/y 



18^(1 - 2A)[(A + 1)Pt] 2 - 2 _ 6A (2 + 3A-9A 2 )(A + l)P r al _ 3A 



Z 4 P?A 



2Z 2 P Q A 



d xyzip + 



2l 2 P 2 \ 

o i ~ 7 

d xj 2 vip 
2~a 2 V 



l 4 Pa 



n = tt - 



67t7[(A + 1)P t ]^ 3 (1 _ 3A) _ (l + 3A)7(A + l)Pr 5 , (1 _ 3A) 
Z 2 P 2 VA ' 2%/aW 



6Vy[(A + l)Pr]^ 3 (1 _ 3A) _ (1 + 3A) v /(A + 1)Pt 5 i (1 _3 A) 

z 2 p 2 Va ' 2VxVv 



12V[(X + 1)P T ] 2 _ 3l _ 2A , (1 + 3A)(A + l)Pr, 2 _3A 



^ 4 P 3 A 



? 2 P a A 



7 2 w 



i 7 4a 3 V i 7i 



d xj 2 ip 



d 3 x^ 2 iptp' 



(B2) 



a = a H 



4aU 



d 27 2 £ 



-2 



Pa = Pa - ^ /d 3 -T^ + ^ / d 3 7 H 2 
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1 _ I 2 P a i _ The intermediary hamiltonian between Ti and 

a 2V reads, 



av. 



(B3) 



Ho 



l 2 P 2 



Pr 



1 



4aV a 3X 2a 



A 



-^t + tt d A x— + — / cParf*v>\i 



2a 



(1 - 3X)PP T - (2+3A ), 



41/ 



+ - 



VA3IV(A + 1)P t 



a -*(5+3A) w + 



3(A+1)P T 



9(A+l)ifr„i- 3 A.,* 



/ 2 P? 



18[(A + l)Pr]^ Q , a _ 9A) _ 3(l_3A)y(ATT)^ a -i(i +3 A) 
Z 2 P Q 3 VA ' 2VVVXP a 



2a 2 V 
' l 4 P 2 



9(1 + 2A)[(A + 1)P T ] _ (1+3A) (2-9A + 9A 2 )/ 2 
2P„ 2 8Fa 2 



ry2y 



J2 1p* 

(B4) 



r 
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